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vasa Background 


¢ During orbital maneuvers and proximity operations, a spacecraft 
fires its thrusters inducing plume impingement loads, heating 
and contamination to itself and to any other nearby spacecraft 

¢ These thruster firings are generally modeled using a combination 
of Computational Fluid Dynamics (CFD) and DSMC simulations 

¢ The Shuttle Plume Impingement Flight EXperiment (SPIFEX) 
produced data that can be compared to a high fidelity simulation 


— Due to the size of the Shuttle thrusters this problem was too resource 
intensive to be solved with DSMC when the experiment flew in 1994 


¢ Objective: 
— Run DSMC Analysis Code (DAC) simulations of specific SPIFEX flight 
test data points 
— Compare the DAC pressure and heating data to the SPIFEX test data 
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Shuttle RCS Thrusters 


CONTROL SYSTEMMODULE primary Thruster (14) 


38 Primary Thrusters (14 Forward, 12 per Aft Pod) 
Daconnect Thrust = 870 Ibs 


Bracket 


6 Vernier Thrusters (2 Forward, 2 per Aft Pod) 
Thrust = 24 Ibs 
"| Propellants: N,O, and MMH 


AFT OMS/RCS POD 


Helium Tanks 


Propellant 
Manifold Valves 


cidizer Tank 
Fual oe | 


Purge and 
Checkout Pane! 


Primary Thrusters 
{12 per Aft Pod) 


https://www.slideshare.net/a 
ticourses/fundamentals-of- 


Reaction Control System A TI rockets-missiles 


Pressurization Components 


Oxidizer Tank 
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Shuttle Plume Impingement 
Flight Experiment (SPIFEX} 


* Flew on STS-64 (Sep. 1994) i 

e Plate with pressure and heat rate 
sensors attached at the end of a 
boom moved around by the 


Shuttle robotic arm al S = 
o> Hard Dia 
¢ Around hundred test points using | ac Hard Disk 


: : é \ Xo 5 
multiple thruster combinations (\ = RCS Hand Controller 


PICS1: Sentran 


¢ Simulated test points: 


— F3U -> Cover range of locations ih tle hnetilal 
in the plume (Centerline versus Pics? —_picse 
high angles) 


— Norm-Z (F3U+L1U+R1U) -> Look 
at interaction zone between 
plumes 


PICS3 PICS! pics2_ PICS11 
/ @ 


4 PICS10 


PICS9 


LMS Plate 


PICS Plate 
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Vacuum Thruster Plume 
Simulations 


¢ Span several flow regimes from continuum inside the nozzle to transitional 
inside the near-field plume to free molecular in the far-field plume 
¢ Currently cannot be modeled with a single solver but must instead use a 
multi-step approach: 
— Use a Computational Fluid Dynamics (CFD) solver in the continuum regions (GASP) 


— Interpolate the CFD solution at some interface to be used as input to the DSMC 
solution 


— UseaDSMC code in the rarefied regions (DAC) 
e Use the Bird breakdown parameter to guide the interface location: 
_ Vx 
B =2— 
cr 
— Interface located where the continuum assumption is valid 


— However, near edge of continuum validity such that DSMC simulation is not too 
computationally expensive 
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SPIFEX Simulations 


e NASA’s DSMC Analysis Code (DAC) 


— Created to solve low density flows such as high altitude plume impingement flows and hypersonic 
reentry flows 


— Parallel, three dimensional code 

— 3D domain meshed using a 2-Level Cartesian grid 

— Usea multi step approach to resolve the flow field 

— Bodies represented using water tight triangulated surfaces 
— Written primarily in FORTRAN with small amount of C 


— Uses the Message Passing Interface (MPI) message passing scheme to effect communication 
between the processors 


¢ SPIFEX simulations parameters: 
— Use nearest neighbor collisions 
— Target a cell size of 2 hard sphere mean free path 
— Target 10 molecules per cell 
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DAC Input Conditions 


¢ Nozzle and near field plume solved with 
the General Aerodynamic Simulation 
Program (GASP) CFD code 
— Chamber Pressure = 152 psi 
— 11 species (CO2, H20, N2, H2, 02, NO, CO, OH, 
N, O, H) and 86 reactions 
¢ Use a Bird breakdown parameter value 


of 0.03 to guide the interface location Nd (1/m") 


1E+23 
¢ Surface data is scaled to match nominal ~ iE.21 
mass flow rate and thrust (see Backup) a 
e Assume a single species with a 
molecular mass of 23.172 g/mol 
(centerline value in the CFD solution) 
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Run Matrix 


Test Case Objective Test33 ~ Test 20 


% 
Normal impingement Test 15 
along plume centerline 


Impingement at Test 81 
intermediate angle of _ - 
attack along plume Test 80 
centerline 


— 


Normal impingement 
at high angle off 
centerline 


Norm-Z Impingement near dual \ 
interaction region 
81 | Impingement near 
triple interaction ~~ 
region > 4 
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Run Statistics 


Test Case Statistics 


5.8B molecules 
359M cells 


5.8B molecules 
356M cells 


Mean free path (m) 


1.00E+00 
1.00E-01 
1.00E-02 


1.00E-03 


6.1B molecules 
379M cells 
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F3U Test 15 


X=8.25m 


P (N/m?) 


Y=0.0m Nd (1/m?) 

2E+22 
2E+21 
2E+20 
2E+19 
2E+18 
Y 2E+17 


Z=28.8m 
~ ~ a; 
we» 


pt 
o ah 
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F3U Test 20 


X=8.25m 


| 
| 
P (N/m?) 
| 
| 


10 


9 
8 
7 
6 
5 


Y=0.0m Nd (1/m?) 

2E+22 
2E+21 
2E+20 
2E+19 
2E+18 
Y 2E+17 


Z=28.8m 
~ ~ a; 
we» 


> > ; 
| 
)) iW a 
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F3U Test 33 


SSNS SSS ——————————————————————————————————————————————————————————————————— 
an Nd (1/m’) 
a 2E+22 
2E+21 
ee ee 2E+20 
2E+19 
2E+18 
: 2E+17 


Z=19.7m 


 ) 


On 
, 
. a Ke ) 
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Nd (1/m’) 

2E+22 
2E+21 
2E+20 
2E+19 
2E+18 
2E+17 
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nes “¢ 


xs 


+ ers Se 
sf. % ; 


e 


P (N/m?) 


Norm-Z Test 81 


P (N/m?) 


/ Nd (1/m’) 
2E+22 


- 2E+21 
2E+20 
2E+19 

— 2E+18 
2E+17 
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Pressure Comparisons 


Good agreement for all cases 

Very good agreement for plume perpendicular to 
the flow field (Test 15) 

ia agreement for plate at an angle of attack 


> @ SPIFEX CAP 
= 8 CDAC CAP 
= @ SPIFEX SEN 
= 6 
a O DAC SEN 
"nn 
= A SPIFEX Fz/A 
= A 

A DAC Fz/A 


Test 15 Test 20 Test 33 Test 80 Test 81 
Test Number 
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Pressure Comparisons 


DAC pressures are within 50% of SPIFEX Capacitance 

Manometer results 

¢ Large uncertainties in SPIFEX solutions for Tests 33 and 
80 

* Note: Tests 33 and 81 pressures are less than 0.5 N/m? 


2.2 
2.0 
1.8 


= 1.6 
1S) 
<x 1.4 
Le CO DAC CAP 
ow 1.2 @ SPIFEX SEN 
v 
s 1.0 © DAC SEN 
nn 
d 0.8 A SPIFEX Fz/A 
A. 
oO 0.6 A DAC Fz/A 
Z 
3 0.4 
~ 02 
0.0 


Test 15 Test 20 Test 33 Test 80 Test 81 
Test Number 
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Forces and Moments 
on LMS Plate 


0.5 
Wm SPIFEX Fx 
e 
> ripe ey Good agreement for 
= largest force 
Y @ SPIFEX Fy 
2 components for 
S O DAC Fy 
normal 
A SPIFEX Fz ae 
impingement 
adel ¢ Worst agreement 
for plate at an angle 
ee of attack (Test 20) 
O Xx 
=E and for Norm-Z 
DAC M 
2 ikiaehs cases (Tests 80 and 
= @ SPIFEX M 
2 Y 81) 
E O DAC My 
oO 
2 A SPIFEX Mz 
A DAC Mz 
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Heat Rates 


¢ Very good comparison for high heat rate cases 
¢ DAC underestimates the heat rates for the high angle off centerline case (Test 33) 


“ “@-Test 15 
7 ~_\AL ™ 
ene A, vee aa 
0.6 ey =@Test 33 


Me a Lines: SPIFEX ——— maeeot 20 
Open Symbols: DAC =@Test 81 


Heat Rate (W/cm2) 
je) 
IN 


0.3 1) Test 15 
— A Test 20 
0.2 — ee — 
¢ ee EE © Test 33 
0.1 — ——__ ~~ N SS 
Soe x Test 80 
0 O © © Q Oo O Q 
O Test 81 
q5 q6 q7 qs q9 qid qil 


Heat Rate Sensor ID 
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Summary 


e F3U: 


— Good agreement for pressure, forces and moments for near normal 
impingement 


— Larger discrepancies for plate at angle of attack and high angle off 
centerline cases 


— Very good agreement for high heat rate cases 
° Norm-Z: 
— Larger discrepancies between DAC and SPIFEX results 
¢ Forward work: 
Rerun the CFD simulation of the nozzle and plume near field 


Add the OMS pods to the shuttle geometry being modeled in the 
Norm-Z simulations 


Run multi species case in DAC 
Run additional test data points 
— Do sensitivity study of the impact of changes in impingement angle 
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Backup 


Scaling to Match Nominal 
Thrust and Mass Flow Rate 


¢ Nominal Values: 
— Mass Flow Rate: 3.1 lbm/s 
— Thrust: 870 lbf 


¢ Final values match the nominal values within 0.5% 


Number Density 21.7% 17.2% 
Scaling 
Velocity -7.7% -6.6% 


Magnitude Scaling 
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Run Statistics 


Test Case Statistics 
5.8B molecules 
359M cells 
5.8B molecules 
356M cells ae 
1.00E+00 
7.94E-01 
6.31E-01 
5.01E-01 
6.1B molecules pees 
2.51E-01 
379M cells SEs: 
1.58E-01 
1.26E-01 
Norm-Z 23.4B molecules 1.00E-01 
1.59B cells 
81 | 22B molecules 
1.4B cells 
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Run Statistics 


Test Case Statistics 


5.8B molecules 
359M cells 


5.8B molecules 
356M cells Population 


6.1B molecules 10 
379M cells 8 


Norm-Z 23.4B molecules 
1.59B cells 
81 | 22B molecules 
1.4B cells 
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